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In our daily lives, after shaking a salad dressing, we see the coarsening of oil droplets 
suspended in vinegar. Such a demixing process is observed everywhere in nature and 
also of technological importance. For a case of high droplet density, domain coars¬ 
ening proceeds with interdroplet collisions and the resulting coalescence. This phe¬ 
nomenon has been explained primarily by the so-called Brownian coagulation mecha¬ 
nism: stochastic thermal forces exerted by molecules induce random motion of indi¬ 
vidual droplets, causing accidental collisions and subsequent interface-tension driven 
coalescence. Contrary to this, we demonstrate that the droplet motion is not random, 
but hydrodynamically driven by the composition Marangoni force due to an interfacial 
tension gradient produced in each droplet as a consequence of composition correla¬ 
tion among droplets. This alters our physical understanding of droplet coarsening in 
immiscible liquid mixtures on a fundamental level. 


Phase separation is one of the most fundamental phys¬ 
ical phenomena that lead to pattern formation in various 
types of systems m, including all kinds of classical con¬ 
densed matter, biological systems (both ordinary [ 7 ] and 
active matter 0 ), quantum systems such as liquid ^He- 
^He mixtures [9] and Bose-Einstein condensates unmi], 
nuclear matter [T^] , and cosmology m- In many of these 
systems, momentum conservation plays a crucial role in 
phase separation. The simplest among such systems is 
a classical binary liquid mixture, where material can be 
transported by both diffusion and hydrodynamic flow. 
Thus, the relevant field variables are the composition of 
one of the components (j) and the fluid flow velocity v. 
The process of phase separation is then described by the 
Ginzburg-Landau-type free energy and a set of cou¬ 
pled dynamical equations of these two variables satisfy¬ 
ing the composition and momentum conservation under 
the incompressiblity condition V • v = 0 (see Methods for 
the explicit forms of the equations). This model is widely 
known as model H in the Hohenberg-Halperin classifica¬ 
tion of dynamical critical phenomena and phase separa¬ 
tion kinetics to dynamic universality classes m- Phase- 
separation behaviour in many of the above-mentioned 
systems can basically be described in a similar frame¬ 
work. Thus, revealing the fundamental physical mecha¬ 
nism of domain coarsening in model H is of significant 
importance for the general understanding of spinodal de¬ 
composition in these systems. 

From a technological point of view, since phase sepa¬ 
ration starts at nanoscale and its characteristic size pro¬ 
ceeds to grow indefinitely with time, we can control the 
domain size of materials intentionally by freezing the pro¬ 
cess at a certain time by using liquid-solid transitions 
such as glass transition and crystallization, by chemical 
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reaction, or by emulsifying liquid droplets with surfac¬ 
tants. This strategy is widely used in industrial processes 
of various materials such as liquid mixtures, emulsions, 
polymer blends, metallic alloys, ceramics, cosmetics, and 
foods. Among various types of phase separation, liquid 
phase separation is quite important in materials science 
since a liquid state is most suitable for materials process¬ 
ing. Thus, the understanding of domain coarsening dy¬ 
namics in the demixing process of liquid mixtures is not 
only of fundamental importance but also of crucial tech¬ 
nological importance for controlling domain structures of 
materials. 

The theoretical framework of model H has a very firm 
basis since it replies solely on the composition and mo¬ 
mentum conservation, and there is little doubt on its va¬ 
lidity. This phase-field model is widely used for numer¬ 
ical simulations not only in physics, but also in materi¬ 
als science and chemical engineering communities. How¬ 
ever, the analytical theoretical analysis of model H is not 
an easy task due to its non-locality, non-linearity, com¬ 
plex non-local dynamical coupling between (f and v, and 
stochastic nature due to thermal noise. Nevertheless, the 
concept of dynamical scaling provides us with a power¬ 
ful scaling argument for the self-similar domain growth 
in the late stage of phase separation after the formation 
of a sharp domain interface HHS [15], on the basis of 
clear physical pictures on the elementary process of do¬ 
main coarsening. Thus, the basic mechanisms of domain 
coarsening in immiscible liquid mixtures are now believed 
to be reasonably understood m- Recently, research in¬ 
terests have shifted towards other aspects of phase sep¬ 
aration including aging dynamics during the coarsening 
process m, geometrical features m, and response func¬ 
tions of coarsening systems [T8|, inertia effects usms], 
and viscoelastic effects [lHoj. 

Here we summarize the current understanding of liq¬ 
uid phase separation. We consider a binary liquid mix¬ 
ture, whose components have the same viscosity r]. A 



2 


schematic phase diagram is shown in Fig. [Ti, together 
with typical phase-separation structures (Fig. ^,c). The 
relevant domain coarsening mechanism depends solely 
on the volume fraction of the minority phase For 
small ^ (case A), the evaporation-condensation (Lifshits- 
Slyozov-Wagner (LSW)) mechanism [2TH23] is respon¬ 
sible for droplet coarsening: steady diffusion flux from 
smaller droplets to nearby larger droplets leads to the 
growth of the latter at the expense of the former. In 
this case, the average domain size (R) grows as {R{t)) = 
[{4/27)Dom"/^ 0 , where t is the time, Dq is the diffu¬ 
sion constant of a component molecule (or atom), and 
^ is the correlation length of composition fluctuations 
(or the interface thickness). In this mechanism, trans¬ 
lational diffusion of component molecules (or atoms) is 
the main transport process. Near a symmetric compo¬ 
sition (^ ^ 1/2) (case B), bicontinuous spinodal de¬ 
composition takes place (Fig. Ef)- There, the hydro- 
dynamic coarsening (Siggia’s) mechanism characteristic 
to a bicontinuous pattern leads to rapid coarsening m- 
{R{t)) ^ 0.1{a/r])t 0, where a is the interface tension 
between the two phases. In this mechanism, hydrody¬ 
namic flow is the main transport process. For interme¬ 
diate ^ (case C), droplet spinodal decomposition takes 
place (Fig. [^). It has been believed that the Brownian- 
coagulation (Binder-Stauffer-Siggia (BSS)) mechanism is 
mainly responsible for the coarsening laiES]: droplets 
grow by accidental collisions between droplets undergo¬ 
ing random Brownian motion, which leads to the coars¬ 
ening law {R{t)) = [{^ksT^/ [H [26], where 
/cbT is the thermal energy. In this mechanism, thermal 
Brownian motion of droplets (rather than molecules) is 
the main transport process. The growth exponents of 
these coarsening laws have been confirmed by both ex¬ 
periments [T7H33] and numerical simulations [191 134H36) . 

The coarsening mechanisms for case A and B are rather 
well established, whereas that for case C is not so clear 
due to difficulties associated with many-body effects and 
hydrodynamic effects. The simplifications of case A and 
B come from the fact that for case A transport is domi¬ 
nated by diffusion and for case B by hydrodynamic flow. 
On the other hand, both may play an important role for 
case C in a complex manner. Although there has been a 
wide belief that the primary coarsening mechanism is the 
BSS mechanism, additional mechanisms have also been 
suggested. For example, Tanaka observed the process 
of droplet coarsening in a binary liquid mixtures con¬ 
fined between two glass plates with optical microscopy, 
and found 1261137 ] that droplet collision enhances sub¬ 
sequent collisions via hydrodynamic flow (phenomenon 
1) and droplets sometimes move directionally, contrary 
to the expectation from random Brownian motion (phe¬ 
nomenon 2). Phenomenon 1 is enhanced for a quasi-two- 
dimensional (quasi-2D) situation due to stronger hydro- 
dynamic interactions [38] , but even for 3D its importance 
has been pointed out when the droplet volume fraction is 
very high [36l|39]. On phenomenon 2, some theories were 
proposed to explain directional motion of droplets by a 



FIG. 1: A schematic phase diagram of a critical bi¬ 
nary mixture with typical phase-separation struc¬ 
tures. a, Schematic phase diagram. In the (light yellow- 
green) metastable region (ie., the low $ region ($ 
0.21)), nucleation-growth-type phase separation takes place 
and there droplets coarsen by the evaporation-condensation 
(or, LSW) mechanism: diffusional transport from smaller to 
larger droplets. Near the symmetric composition (the light 
blue region, or, $ ~ 1/2), bicontinuous spinodal decomposi¬ 
tion (b) takes place and the domain coarsening is governed 
by the Siggia’s hydrodynamic coarsening mechanism: hydro- 
dynamic transport from thin to thick parts of bicontinuous 
tube structures. In the intermediate composition region (the 
green region), droplet spinodal decomposition (c) takes place 
and its coarsening mechanism has been believed to be the 
Brownian coagulation mechanism: random Brownian motion 
of droplets and the resulting collisions. Contrary to this, we 
propose a new coarsening mechanism, in which droplet mo¬ 
tion is induced by the composition Marangoni effect due to 
the diffusional composition coupling among droplets, and not 
by random Brownian motion. The border between bicontin¬ 
uous and droplet spinodal decomposition is located around 
$ = 0.34, according to our simulation. 


coupling between composition and velocity fields [40l|4T]; 
however, they turned out to be wrong even on a quali¬ 
tative level, as will be described later. The difficulty of 
this problem stems from the non-local nature of diffusion 
and the many-body and long-range nature of its coupling 
with hydrodynamic degrees of freedom. Although direc¬ 
tional motion of droplets was experimentally reported 
for such a quasi-two-dimensional situation |26l |37j , there 
has so far been no firm experimental evidence supporting 
such spontaneous motion of droplets for a bulk 3D sys¬ 
tem. Because of this lack of clear experimental evidence 
together with the theoretical difficulty mentioned above, 
phenomenon (ii) has not attracted much attention. Thus, 
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it is now widely believed that droplet coarsening in im¬ 
miscible liquid mixtures can be explained by the BSS 
mechanism HHS]. 

In this Article, we will demonstrate by numerical simu¬ 
lations that, contrary to the BSS mechanism, droplet mo¬ 
tion is not random, but hydrodynamically driven by the 
composition Marangoni force, which is induced by long- 
range diffusional composition correlation among droplets, 
indicating that a dynamical coupling between the com¬ 
position and velocity field leads to a novel mechanism of 
efficient droplet coarsening in droplet spinodal decompo¬ 
sition. 

Results 

Phase diagram and the setting of basic parame¬ 
ters 

Before presenting results of phase-separation simula¬ 
tion, we summarize phase separation behaviour of a bi¬ 
nary fluid mixture, focusing on how it depends on the 
composition (j) and the temperature T (see the schematic 
phase diagram. Fig. Ef)- Here we consider phase demix¬ 
ing induced by an instantaneous temperature quench at 
time t = 0 from an initial temperature in the one-phase 
region (white region) to a final temperature T in the two- 
phase region. First the phase separation behaviour can 
be grouped into two types: nucleation&growth (NG)- 
type phase separation in the metastable region (the light 
yellow region in Fig. Ef) and spinodal-decomposition 
(SD)-type one in the unstable region. Then, SD-type 
phase separation can further be classified into droplet and 
bicontinuous SD. Domain coarsening is mainly driven 
by the LSW mechanism for NG-type phase separation, 
which takes place for ^ ^ 1, and by Siggia’s mecha¬ 
nism for bicontinuous SD, which takes place for $ ^ 1/2. 
For droplet SD, which takes place for the intermediate ^ 
Houghly, r\j 0.21 < ^ 0.34; the green region in Fig. 

[^), it has been widely believed that the BSS mechanism 
is responsible for domain coarsening. We note that the 
boundary between NG- and SD-type phase separation is 
sharp only in the mean-field limit, and should be diffuse 
for ordinary binary liquid mixtures [42]. 

The composition region on which we focus here is the 
region of droplet SD. The thermodynamic state of a bi¬ 
nary fluid mixture is characterized by the composition 
(j) and the temperature difference between the critical 
temperature Tc and the temperature T. Then the cor¬ 
relation length of the composition fluctuation which 
also characterizes the interface thickness, is described 
as ^ = a|(T — Tc)/Tc|“^, where the critical exponent 
z/ ^ 0.63 for 3D Ising universality class, to which a binary 
mixture belongs HE]. On the other hand, the kinetics 
is characterized by the fluid viscosity r]. Then the cooper¬ 
ative diffusion constant is expressed as = kBT/{6TTriC) 
and the characteristic lifetime of composition fluctuations 
is given by /ksT. As shown in Meth¬ 

ods, the phase separation behaviour of a binary fluid mix¬ 
ture can be described solely by four dimensionless param¬ 
eters. One is the scaled composition (j) = where 


is the final equilibrium composition, or equivalently the 
volume fraction of the minority phase 4>. The others are 
the fluidity parameter M, which is a measure of the rel¬ 
ative importance of the hydrodynamic to the diffusional 
transport [43], the thermal noise parameter S, which is 
the strength of the noise term (see Methods), and the 
Reynolds number Re, which characterizes the relative im¬ 
portance of the inertia term to the viscous term. A and 
B are defined respectively as A = 67r|7|^^0g = 

and B = [35] [43] . We note that a in 

these relations is the interface tension in an equilibrium 
two-liquid coexistence state (see Methods for its expres¬ 
sion in the Ginzburg-Landau model). We set Re = 0 
since in ordinary droplet phase separation we can safely 
neglect the inertia effects m (see Methods). 

We note that near a critical point Tc, it is known 
that the renormalized value of M is a universal constant, 
M = 4 ^ 8 [35j|43], since the equilibrium interface tension 
is given by a = (0.2 ^ 0.4 )/cbT/^^ according to the two- 
scale-factor universality [4|. However, far from a critical 
point, this parameter can become very large for typical 
immiscible binary liquids (e.^., for a water/hexane mix¬ 
ture, A ^ 100), whereas it becomes very small for poly¬ 
mer mixtures since A ^ N~^ {N: the degree of poly¬ 
merization of polymer) [H [35]. On the other hand, the 
estimation of the value of B needs some care, since this 
represents the strength of thermal fluctuations that di¬ 
rectly lead to the break-down of the mean-field picture 
due to renormalization effects: A larger value of B shifts 
Tc to a lower temperature. We estimate that B should 
be the order of 1 near a critical point, and it decreases 
when the temperature is much lower than Tc. So, in 
our simulations, we chose A = b and S = 1 as typical 
parameters for a system near a critical point. We note 
that the choice of these values of A and B does not af¬ 
fect our basic conclusion on a qualitative level. As shown 
below, our new mechanism, where hydrodynamic trans¬ 
port plays a key role, becomes more important with an 
increase in the distance from a critical point, ie., for 
larger A and smaller B. Note that A/B is a measure 
of the ratio of the hydrodynamic transport to the noise- 
driven diffusional transport. This means that our new 
mechanism should become more important for a liquid 
mixture far from a critical point, which is often the case 
for realistic applications, than for a critical mixture. 

Finally, we note that our study concerns only the so- 
called late stage of phase separation after the formation 
of a sharp domain interface, where the domain size is 
large enough compared to the interface thickness: 

Below we study 4> = 0.25, which is in the composition 
region of droplet SD (see Fig. Ef)- Hereafter we use the 
scaled domain size R = R/^ and the scaled time i = tjr^ 
to describe the coarsening dynamics. 

Results without thermal noise 

Firstly, we show a coarsening process during droplet 
phase separation in Fig. ii, which is simulated by solving 
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the dynamical equations of model H (see Methods) with¬ 
out thermal noise (see Methods). Because of the absence 
of noise (z.e., B = 0), the droplet interface is smooth. 
The time evolution of the average size {R) and the total 
number of droplets N for the fluidity parameter A = b 
(black) and M = 50 (red) are shown in Fig. and c, re¬ 
spectively. In addition to M = 5, we use the large fluidity 
parameter M = 50 to access a large separation between 
the interfacial thickness and the droplet size as well as to 
see hydrodynamics effects clearly. Because of the com¬ 
position asymmetry, the symmetry of initial composition 
fluctuations around the average composition, which holds 
in the initial linear regime, is broken once the non-linear 
terms in the free energy start to come into play, and thus 
many small droplets are formed by fragmentation in the 
early stage. As stated above, in this study we are con¬ 
cerned only with the late stage of coarsening after the 
formation of a sharp domain interface. In the late stage, 
we find (R) oc and N oc for both fluidity param¬ 
eters (Fig. |^,c). Even when thermal noise is absent, we 
can see spontaneous motion of droplets and subsequent 
collisions: the droplet coarsening proceeds even without 
thermal noise. This result can never be explained by 
the Brownian coagulation mechanism, in which thermal 
noise is the only cause of droplet motion and coarsening. 
This is the most direct evidence for the presence of an un¬ 
known coarsening mechanism other than the BSS mech¬ 
anism, which is responsible for the spontaneous motion 
of droplets and their collisions. Because of the absence of 
thermal noise, the droplet coarsening process is perfectly 
deterministic in this case, and thus the initial compo¬ 
sition noise introduced at t = 0 completely determines 
what takes place afterwards. 

For M = 5, we can see both collisions of droplets and 
gradual disappearance of droplets due to evaporation, al¬ 
though the former is much more frequent than the latter. 
For A = 50, on the other hand, the growth of droplets 
proceeds almost solely through direct collisions between 
droplets and evaporation of droplets is very rare. This 
difference can be explained by the fact that M is a mea¬ 
sure of the relative importance of hydrodynamic to dif- 
fusional transport. Figures and f show respectively 
the average diffusion flux and the average magnitude of 
the velocity of droplets for A = bO during a coarsening 
process. Each spike reflects interdroplet collision. We 
can see that, in the late stage, collisions take place af¬ 
ter a complete decay of flow and thus are not affected 
by the preceding collision even for M = 50, clearly indi¬ 
cating that the migration of droplets and the following 
collision are due to spontaneous motions of droplets and 
‘not’ induced directly by the preceding collision. 

Results with thermal noise 

Next we show a coarsening process during droplet 
phase separation with thermal noise in both composi¬ 
tion and force fields. We note that the incorporation of 
both types of noise in the model H equations has been 
technically challenging j44ll46] . However, this has been 


properly done for 2D in the framework of model H by 
Camley et al. [44]. We have employed their method to 
simulate a 3D system (see Methods). The merit of this 
method is that we can drop the inertia term. 

The process of phase demixing is shown in Eig. (see 
also Supplementary Movie 1). The simulation is made for 
the fluidity parameter A = b and the thermal noise pa¬ 
rameter B = 1 (see the above for the choices of these val¬ 
ues). We can clearly see the temporal growth of droplets, 
which are induced by spontaneous droplet motion and 
the resulting collision and coalescence, as observed in Eig. 
§1. Unlike in Eig. however, the interface of droplets 
are rough due to thermal fluctuation effects. 

As shown in Eig. [^, the structure factor S{k) can 
be scaled by using the characteristic wavenumber ki 
(see Methods), indicating that droplet patterns grow 
self-similarly The time evolutions of the average 

droplet size (R) and the number of droplets N are shown 
in Eig. [^. We find (R) ^ and N ^ The ex¬ 
ponents are consistent with the prediction of the BSS 
mechanism, where the coarsening rate is predicted as 
{R{t))^ = {b^kBT/b7rr])t jH |26]. In our scaled units, 
this relation becomes {R{i))^ = ((34>MS)/57r)t = 0.24t 
for M = 5 and B = 1. However, what we find is 
{R{i))^ = 1.6t. This indicates that the droplet size (R) 
grows 1.9 times faster than the prediction of the BSS 
mechanism although the exponent of the power law is 
the same. We note that our simulation result is much 
consistent with the experimental finding of Perrot et al. 
under microgravity [31] that the droplet size grows twice 
faster than the prediction of the BSS mechanism. This 
supports the validity of our simulations with the above 
choices of A and B. 

Spontaneous directional motion of individual 
droplets 

To seek the cause of the non-Brownian nature of 
droplet coarsening in simulations without noise (Eig. 
and the faster coarsening than the BSS prediction in sim¬ 
ulations with noise (Eig. |^, we analyse trajectories of 
droplets during the interval of two successive collisions 
for the case with thermal noise. Eirst we focus on the 
nature of droplet motion, z.e., whether it is completely 
random as assumed in the BSS mechanism (Eig. or 
directional (Eig. Ef)- We carefully choose the observa¬ 
tion period to avoid the effects of flow induced by colli¬ 
sions of droplets and to see only spontaneous motion of 
the droplets. In Eig. H, a typical trajectory of a droplet 
during a certain time period are shown. We can clearly 
see that the droplet moves directionally (towards a neigh¬ 
bouring droplet). We confirm that this type of motion 
is common to most of droplets and not special (see, e.g.. 
Fig. |§i). Below we will show that it is this directional 
motion that leads to the fast coarsening process. We note 
that this motion is not due to inertial effects [36] since 
our simulations are performed under the Stokes approx¬ 
imation {Re = 0, or without the inertia term), which is 
valid for ordinary droplet phase separation. Thus, there 
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FIG. 2: Coarsening dynamics of droplet phase separation without thermal noise, a, Coarsening process of droplet 
spinodal decomposition at $ = 0.25 for ^ = 5 and 8 = 0. Blue and red coloured surfaces show front and back sides of a droplet 
interface, respectively. Time evolution of the average droplet size {R{t)) and the total number of droplet N(t) are shown in 
b and c for ^ = 5 (black) and ^ = 50 (red), respectively. The green straight lines have the slope of 1/3 for a and -1 for b. 
d, Temporal change in the average diffusion flux (j) for A = 5. Note that the characteristic decay time of the diffusion flux 
after a collision is ~ O'KrjR^a/k^T. In our unit, this time scale is roughly proportional to W (e.g., ~ 400 for R = 20). e, 
Temporal change in the average velocity magnitude (y) for ^ = 5. Note that the characteristic decay time of the velocity field 
after a collision is Tv ~ riR/a, which is the characteristic shape recovery time of a deformed droplet to a spherical shape. In our 
unit, this time scale is roughly proportional to R/A (e.g., ~ 4 for ^ = 20). Each sharp spike in d and e reflects interdroplet 
collision. The insets in d and e are the same plots but with a semi-log x axis. 



FIG. 3: Coarsening dynamics of droplet phase separation with thermal noise, a, Temporal change in the droplet 
structure in the late stages of the phase-separation process at 4> = 0.25 for A = 5 and 8=1. Blue and red coloured surfaces 
show front and back sides of a droplet interface, respectively, b, Scaling plots of the normalised structure factor S{k). The 
collapse onto the master curve indicates the self-similar nature of pattern evolution, c, Time evolution of the average droplet 
size {R(t)) and the total number of droplets N(t) during phase separation. The blue and green lines have the slope of -1 and 
1/3, respectively. 


must be some thermodynamic force acting on droplets For a droplet of radius R which is undergoing ran- 
that drives droplets to spontaneously move directionally. dom Brownian motion freely, its mean-square displace- 
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FIG. 4: Evidence for directional motion of droplets. 

Schematic pictures for random and directional motion are 
shown in a and b, respectively. Case a is for non-interacting 
droplets and case b for interacting droplets, c: An example 
of a trajectory of the centre of mass of a single droplet during 
i = 12120 — 13520 in the process shown in Fig. 3. The arrow 
indicates the direction of the motion. The result clearly shows 
directional translational motion in the phase-separating pro¬ 
cess, strongly suggesting that the droplet motion is the type 
depicted in b and not in a. d: An example of displacements 


normalized by radius, y R(t)5r(t)^ of droplets during a time 

interval (open symbols: i = 11320 — 11720; filled symbols: 
t = 12120 — 13520) in the process shown in Fig. 3. We note 
that during this time interval no collisions take place. Sym¬ 
bols with different colours are for different droplets in the 
simulation box. The blue straight line is \/R6r = 0.4?^^ (see 
Eq. (1)) and the red dashed line has the slope of 1. The 
displacements are always much larger than the prediction of 
Brownian motion (the blue straight line) for non-evaporating 
droplets. Note that droplets having displacements smaller 
than the prediction are evaporating ones. 


ment, (Ar^), should be given by (Ar(t)^) = = 

{ksT/{b7rr]R))t. In our scaled units, for ^ = 5 and 

B = h 


^ 0.4f/2. 

IOtt 


( 1 ) 


In order to compare the displacements of droplets with 

different radii, we plot ^/rSt for all droplets in Fig. |^. 
In this plot, the difference in the droplet size are scaled 
out; and, thus, the displacements of all droplets should 
agree to the prediction of Eq. 0 and thus fall on the 
blue solid line in Fig. ii, if droplets are doing free ran¬ 
dom Brownian motion. However, we find that most of 
droplets have much larger displacements far beyond the 
prediction of Eq. 0 (or, the blue solid lines), strongly 
indicating the non-random directional nature of droplet 


motion. During this time interval of observation, the 
change in the droplet radius is within a few %, so this 
effect does not play any role. Hydrodynamic interactions 
between droplets and finite-size effects should slow down 
the self-diffusion of droplets. Furthermore, the displace¬ 
ment should still increase in proportion to \/t as long 
as the motion is random. Nonetheless, Fig. H clearly 
show that the displacements of droplets are much larger 
than the prediction of Eq. Q and even increase almost 
linearly at later times. This cannot be explained by the 
scenario based on random Brownian motion. 

Without noise, we should be able to see clearly the 
relation between the thermodynamic force and the re¬ 
sulting transport, diffusional and hydrodynamic, with¬ 
out suffering from noise that obscures the details. Thus, 
we study the chemical potential on droplet interfaces for 
the case without thermal noise, and the result at a cer¬ 
tain time is shown in Fig. 0- We can clearly see that the 
chemical potential on the interface of a droplet, or the in¬ 
terfacial composition profile, is not spherical symmetric 
but rather anisotropic, reflecting the spatial configura¬ 
tion of the neighbouring droplets. As will be shown be¬ 
low, this causes the composition Marangoni force, which 
induces spontaneous directional motion of droplets. Nev¬ 
ertheless, the droplets are nearly spherical because of 
the action of the Laplace pressure Ap = 2(t j R , which 
induces a quick hydrodynamic relaxation to the stable 
spherical geometry with a characteristic relaxation time 
of Ty ^ r]R/a. A small variation of the interface tension. 
Act, in a droplet does not cause shape deformation of the 
droplet, but causes a drastic change in the type of droplet 
motion from stochastic to deterministic nature. We note 
that the former is a relative effect, whereas the latter is 
an absolute effect. 

Figure is an enlargement of the blue small droplet 
and its surrounding droplets shown in Fig. 0. The ar¬ 
rows show the flow field at their interfaces. The flow field 
is caused by an imbalance between the Laplace pressure 
and the thermodynamic force under the incompressibil¬ 
ity condition: the break down of the spherical symmetry 
of the interfacial chemical potential in a droplet and the 
resulting tangential force acting along the interface allow 
the directional flow field to be induced even under the 
condition of V-v = 0, unlike the case of an isolated spher¬ 
ical droplet. We can see that the smaller bluish droplet in 
a higher free energy state is slowly evaporating and the 
neighbouring droplets in lower energy states are grow¬ 
ing. Thus the interface of a neighbouring droplet that is 
closer to the evaporating small droplet has smaller chem¬ 
ical potential (more bluish) than the opposite side. This 
anisotropy of the chemical potential on the interface of 
each droplet is the origin of the deterministic directional 
motion of the droplet. The variation of the chemical po¬ 
tential on a droplet interface is less than a few %, but the 
resulting composition Marangoni force is strong enough 
to overwhelm the thermal force acting on a droplet. 

The emerging physical picture is as follows. An 
small (evaporating) droplet causes deterministic flow 
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FIG. 5: Coupling between composition and velocity 
fields, a, A snapshot of droplets during phase separation 
without thermal noise. The colour reflects differences of 
the chemical potential on the interface. Reddish and bluish 
colour means higher and lower chemical potential, respec¬ 
tively. As shown in the colour bar, the colour is assigned 
for each droplet, which is bound by the minimum and max¬ 
imum value of the chemical potential /x. This is because the 
amplitude of the change of /x inside each droplet interface is 
much smaller than the fluctuations of the absolute value of 
/X of each droplet. We note that the magnitude of difference 
relative to mean value is less than 3% even for a droplet which 
has the largest intradroplet variation, b, Enlargement of the 
blue small droplet in a and its surrounding droplets. The 
meaning of colour is the same as in a. The velocity field at 
the droplet interfaces is also shown by the arrows. A longer 
arrow means high velocity. 


fields on its neighbouring larger droplets, which lead 
to motion of the neighbouring droplets towards it 
and eventually induce an interdroplet collision, if the 
droplet does not evaporate before collision. In other 
words, a larger droplet tends to eat a neighbouring 
smaller droplet by direct collision. This hydrodynamic 
eating speed is much faster than the diffusional one 
in the evaporation-condensation (LSW) mechanism, as 
long as hydrodynamic transport is faster than diffu¬ 
sional transport. More precisely, however, the com¬ 
position correlation between neighbouring droplets in¬ 
duces both the hydrodynamic and diffusional transport. 
Thus, the crossover between the present mechanism and 
the evaporation-condensation mechanism should be con¬ 
trolled by the competition between hydrodynamic and 
diffusional transport. The former is always dominant as 
long as is not so low. For lower however, the di¬ 
rect compositional correlation between droplets, which is 
the origin of the thermodynamic force inducing hydro- 
dynamic droplet motion, becomes weaker due to the less 
overlap of the composition fields around droplets. For 
very low thus, droplets interact with each other via the 
mean-field matrix, as assumed in the original mean-field 
theory of the LSW mechanism Enin]- In this regime, 
the LSW mechanism becomes dominant. We speculate 
that the crossover takes place around ^ ~ 0.21 separating 
NG and SD (see Fig. since by definition droplets are 


formed rather independently for NG whereas in a corre¬ 
lated manner for SD. But this crossover composition may 
also depend on the value of A. Although this is a very 
interesting problem, we leave it for future investigation. 

Analysis of a pair of droplets 

To elucidate the exact physical mechanism, we also 
study the motion of pairs of droplets (see Methods). Typ¬ 
ical results are shown in Fig. I§i and b, where the chemi¬ 
cal potential and the flow field on the interfaces are shown 
as in Fig. [^. Figure and d show the composition pro¬ 
file of two droplets along the line connecting the centres 
of mass of the two droplets for case a and b. Figure 
and f are their enlargements, which clearly show the pres¬ 
ence of the composition difference by A0 between the two 
sides of a droplet. We can see that such an asymmetric 
interfacial profile of a droplet induces a chemical poten¬ 
tial gradient in the droplet and leads to spontaneous di¬ 
rectional droplet motion. Finally, the temporal change 
in the centre-of-mass velocities of the two droplets are 
shown in Fig. and h, respectively, for case a and b. 

For the case of Fig. |^, we calculate the interfacial ten¬ 
sion (j for each interface by integrating 
m across the interface along the radial direction. A part 
of this thermodynamic force is balanced by the Laplace 
pressure Ap = 2ajR. For the left droplet whose ra¬ 
dius is 11.5 and velocity is -0.0045, the interfacial ten¬ 
sion is 0.9286 for its left side and 0.9296 for the right 
side. This small variation of the interfacial tension is 
consistent with a nearly spherical shape of droplets, as 
discussed above. It is the interfacial tension gradient due 
to the composition gradient on the interface that leads to 
the motion of a droplet. This is known as the composition 
Marangoni effect imiiH]. We note that it is well known 
that this effect induces spontaneous droplet motion for 
droplets with weak surfactants [49] and those subjected 
to composition or temperature gradients [50l [51]. We 
stress that unlike these well-known cases, for the present 
case the interfacial tension gradient in a droplet is pro¬ 
duced by a non-trivial coupling between the composition 
fields around droplets. The analytical expression for the 
droplet velocity under a uniform composition gradient is 
give as V = —{R/{br])){da/dx) [47]. By inserting to this 
relation the above values of a estimated from the simu¬ 
lation, we estimate the droplet velocity as |V| ^ 0.0050, 
which is much consistent with the directly measured one 
(|V| ^ 0.0045) (see Fig. [^). A slight difference may 
come from the difference in the geometry, z.e., a single 
particle under a uniform composition gradient (theory) 
vs. a droplet pair under diffusional coupling (our sim¬ 
ulation), and the finite size of the simulation box. To 
avoid the finite-size effects, we need to perform larger size 
simulations, but unfortunately a large numerical cost of 
hydrodynamic simulations do not allow us to do so. The 
velocities of the droplets in Fig. 7b are also estimated 
from the same analysis, and they also agree well with 
the measured ones (see Fig. 7b and its caption). These 
results strongly support our mechanism. 
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FIG. 6: Motion of a pair of growing droplets, a, Motion of two growing droplets of the same size. The centre-of-mass 
distance between the two droplets is 35 and the initial radii of the two droplets are both 10. The velocity field and the chemical 
potential at the droplet interfaces are also shown, as in Fig. The interfacial tension of the left droplet is 0.9286 for its left 
side and 0.9296 for the right side around i = 100. The theoretically estimated velocity from these values is -0.0050, which is 
consistent with the measured velocity. The colour bar is the same as in Fig. We can see the higher interface tension at the 
side of a droplet near another one. This leads to the motion of the droplets increasing their separation, due to composition 
Marangoni effects, b, Motion of a pair of droplets with different size. The centre-of-mass distance between the two droplets is 
45, and the initial radii of the left and right droplets are 10 and 15 respectively. The sides of the droplets nearest each other 
have lower interfacial tension and the two droplets approach each other. The interface tension of the left and right side of the 
left droplet is 0.9303 and 0.9294 respectively and the velocity estimated from these values is +0.0045, whereas the interface 
tension of the left and right side of the right droplet is 0.9289 and 0.9313 respectively and the velocity estimated from these is 
-0.012, around t = 100. These velocities estimated around t = 100 are consistent with the measured ones. The simulation box 
size is 192^ for both a and b. Because of its small size, the results suffer from finite-size effects. c,d: The composition profile 
of the droplets along the line connecting their centres of mass for a and b, respectively. e,f: Enlargement of the composition 
profile of the left droplet shown in c and d, respectively. g,h: The temporal change of the centre-of-mass velocity of the two 
droplets for a and b, respectively. The final steep change of the velocity in g is due to a direct overlap of the interfaces of the 
two droplets upon collision. 


In the particular configuration in Fig. [^, due to the 
perfect symmetry around the centre of mass of the two 
droplets and the absence of noise, coarsening does not 
proceed. We note that the same situation should be re¬ 
alized for the LSW mechanism. This is not the case for 
the asymmetric situation in Fig. where two parti¬ 
cles approach each other and eventually collide. For real 
droplet phase separation, the randomness of droplet size 
and many-body interactions always lead to such asym¬ 
metric situations and their continuous presence is respon¬ 
sible for continuous coarsening. We emphasize that the 
randomness is a direct consequence of thermal fluctu¬ 
ations and thus intrinsic to any phase separation phe¬ 
nomena. Combining this result on pairs of droplets and 
the results shown in Fig. and b, we can conclude 
that the composition Marangoni effect stemming from 
the randomness of droplet size is responsible for sponta¬ 
neous directional motion of droplets and the resulting co¬ 
alescence. The importance of randomness indicates that 
there should be a non-trivial feedback between the coars¬ 
ening and the droplet size distribution. We find that the 
normalized droplet size distributions at various phase- 
separation times can be superimposed on the universal 


curve (see Fig. [^, consistent with the self-similar na¬ 
ture of pattern evolution. We speculate that more fre¬ 
quent collisions between particles with large size differ¬ 
ence may be responsible for the self-similarity, but this 
problem should be studied carefully in the future. Our 
mechanism cannot lead to domain coarsening for a mono- 
disperse emulsion. However, the Brownian coagulation 
mechanism should work even in such a situation; and 
once a collision happens and produces the droplet size 
difference, our mechanism should become operative. 

Roles of other non-trivial mechanisms 

Finally we briefly consider other possibilities. In our 
previous papers ESI |37], we suggested two other non¬ 
trivial droplet collision mechanisms: multiple-collision 
and collision-induced collision mechanisms, which was 
confirmed experimentally and numerically [38l |39l [52] . 
Here we show that neither of them can explain the above- 
described phenomena. The former mechanism is a rather 
obvious one. When two droplets collide with each other, 
the shape of the droplet interface is strongly deformed 
and relaxes to a spherical shape, accompanied by the 
interfacial-tension driven flow. This flow induces hydro- 
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FIG. 7: The temporal change in the normalized 
droplet size distribution function. We can see that 
the normalized droplet size distribution functions at various 
phase-separation times can be scaled onto a universal curve 
at least approximately, consistent with the self-similar nature 
of the domain growth (see Fig. 3b). 


dynamic motion of the surrounding droplets, which leads 
to subsequent collisions. The latter mechanism is, on the 
other hand, a result of a coupling between droplet col¬ 
lision and composition diffusion. After droplet collision, 
strong diffusion flux towards the resulting merged droplet 
is created due to what we call the interface quench effect 
[43] (see Fig. Hf), which is the process of local equi¬ 
libration of the droplet with the matrix. This strong 
diffusion affects the translational motion of other neigh¬ 
bouring droplets due to the dynamical coupling of the 
diffusion and flow fields, which leads to subsequent col¬ 
lisions. Both mechanisms are expected to enhance suc¬ 
cessive collisions, but are effective only immediately after 
previous collisions. To reveal whether the coarsening be¬ 
haviour can be described by these mechanisms or not, 
we calculate the average diffusion flux and the average 
magnitude of the velocity during a coarsening process 
without noise respectively in Fig. and e. Occasional 
interdroplet collisions accompany sharp spikes and sub¬ 
sequent relaxation for both quantities. The characteris¬ 
tic decay times of diffusion flux and flow are given by 
T 0 = R^/Dq = 67rr]R‘^a/{kBT) and Ty = r\RjG^ respec¬ 
tively. As can be seen in Fig. Hi and e, in the late stage 
a collision takes place after the diffusional and hydrody¬ 
namic fluxes induced by the preceding collision have com¬ 
pletely decayed, suggesting that it is not affected directly 
by the preceding collision: there are no chain collisions. 
This strongly indicates that the coarsening behaviour we 
observe in the late stage of phase separation can be ex¬ 
plained by neither of the above-mentioned mechanisms. 

Discussion 

The dynamical scaling argument on domain coarsening 
in the late stage of phase separation has been so success¬ 
ful [IHS]. In this argument, it has been assumed that 
the compositions of the two phases already reach the fi¬ 
nal equilibrium one (\)q in the late stage. For droplet 


phase separation, thus, it is implicitly assumed that the 
inter facial tension a is uniform over a droplet’s interface. 
In liquid mixtures, which are incompressible in ordinary 
conditions, the hydrodynamic flow can be induced only 
by the interfacial force, which acts normal to the inter¬ 
face in the above assumption. For a spherical droplet, 
this leads to the conclusion that the spherical symmet¬ 
ric force fields which act along the interface normal unit 
vector n of a droplet cannot cause any flow under the con¬ 
straint of V • V = 0 and accordingly the interfacial force 
should be balanced by the Laplace pressure Ap = 2g jR 
{a being the interfacial tension) [20] . This is the basis 
for the long belief that droplet motion is not induced 
hydrodynamically but solely by its thermal diffusion. 

In the above, we have shown that this physical pic¬ 
ture is not correct and the interfacial tension is actually 
inhomogeneous over each droplet’s surface due to dif¬ 
fusional composition correlation among droplets, which 
allows interfacial-tension-driven flow even for a spheri¬ 
cal droplet. The key is the fact that droplets exchange 
the component molecules (or atoms) with neighbour¬ 
ing droplets via the surrounding matrix phase by dif¬ 
fusion. The time necessary for a component molecule 
(or atom) (size a) to diffuse over the same distance is 
Ts = L‘^/Dq. On the other hand, the average time neces¬ 
sary for a droplet of radius R to diffuse over an average 
inter-droplet distance L is estimated as td = L^IDr, 
where Dr is the diffusion constant of a droplet. Since 
Do and Dr are inversely proportional to the size of a 
component molecule (or atom), a, and the domain size, 
i7, respectively, the ratio td/ts becomes very large in 
the late stage of droplet SD, where ^ a. Thus the 
composition correlation between neighbouring droplets 
develops via molecular (or atomic) diffusion long before 
an accidental collision between them by thermal Brow¬ 
nian motion takes place. It is this composition correla¬ 
tion among neighbouring droplets that is responsible for 
internal inhomogeneity of the interfacial tension over a 
droplet’s interface and the resulting Marangoni force. 

Here we note [24[ |26] that van der Waals attraction 
between droplets is too weak to affect the motion of 
droplets in the late stage. Furthermore, directional mo¬ 
tion of droplets observed in our previous experiments in a 
quasi-2D situation, where neither gravitational force nor 
electrostatic interactions are relevant, can be explained 
by the present mechanism. 

Here it is worth mentioning that the mechanism we 
find here is essentially different from mechanisms pre¬ 
viously proposed by one of the authors (H.T.) |40| and 
Kumaran m- In the above, we reveal that the force 
acting on a droplet comes from the non-uniformity in the 
amount of the composition jump across the interface in 
the droplet and thus is localized on the interface, whose 
characteristic width is given by This non-uniformity 
of the interfacial tension due to the coupling of the com¬ 
position fields around droplets causes the composition 
Marangoni flow, which leads to the spontaneous motion 
of droplets. In both of the previous theories, however. 
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the non-uniformity of the Gibbs-Thomson condition on 
the droplet interface was completely ignored and only 
the composition gradient in the matrix phase was con¬ 
sidered. The characteristic lengthscale of this inhomoge¬ 
neous composition is the order of R. Thus, the strength 
of the composition Marangoni force in the new mech¬ 
anism is much stronger than the composition-gradient- 
induced force in the previous ones since 1/^ ^ 1/R. 
Furthermore, there is another important, even qualita¬ 
tive, difference between them; for example, both of the 
previous theories [40l [41] (wrongly) predict that for a pair 
of large (growing) droplets or for a pair of small (shrink¬ 
ing) droplets, the two droplets should come towards each 
other, whereas for a pair of large (growing) and small 
(shrinking) droplets the two droplets should move away 
from each other. We stress that this is opposite to what 
is shown in Fig. I§i and b. 

The novel mechanism reported here can be explained 
on an intuitive level as follows. There is a difference in the 
local free energy between droplets, reflecting the differ¬ 
ence in the size between them. Larger droplets are more 
stable and in a lower free energy state simply because 
the interface energy cost is smaller. This leads to the dif¬ 
fusion flux from small droplets to nearby larger droplets 
due to the chemical potential difference. This correlation 
of the composition field between droplets evolves much 
faster than any significant diffusional centre-of-mass mo¬ 
tion of the droplets. This is exactly the same process 
as that of the LSW mechanism: if there are neither hy¬ 
drodynamic degrees of freedom nor translational droplet 
diffusion, this leads to the coarsening of the LSW type. 
However, the presence of hydrodynamic degrees of free¬ 
dom in a system entirely changes the kinetic route for 
lowering the free energy. An intradroplet gradient of 
the interfacial tension generates hydrodynamic motion 
of droplets towards the side with lower interfacial ten¬ 
sion. Thus, larger droplets tend to eat nearby small 
droplets, as in the case of the LSW mechanism. However, 
the crucial difference comes from the fact that the trans¬ 
port process is dominantly controlled by hydrodynamic 
translational motion of droplets due to the composition 
Marangoni effect in our mechanism, whereas by trans¬ 
lational diffusion of molecules (or atoms) in the LSW 
mechanism. 

In a more general perspective, we may say that domain 
coarsening is under the influence of the free energy at any 
moment. This has a significant implication for the uni¬ 
fied physical picture of domain coarsening: all coarsen¬ 
ing mechanisms, including the evaporation-condensation 
mechanism (for 4> ^ 1 ), the hydrodynamic coarsening 
mechanism ( for 4> ^ 1 / 2 ), and the mechanism found 
here (for intermediate 4>), have a common principle of 
material transport: Transport can be either diffusion, 
flow, or by their coupling, but always takes place from 
smaller to larger domains: Domain coarsening obeys a 
rule that “big always wins over small”. We note that 
there is no such a rule for the BSS mechanism, in which 
the process is purely random and stochastic. 


Finally, we stress that our new mechanism may play 
an important role in phase separation of many binary 
mixtures in practical use, which takes place in a state 
far from a critical point. We note that our mechanism 
should be more operative for a mixture of lower viscosity, 
as mentioned earlier. We also stress that our mechanism 
can operate even near zero temperature, as shown in Fig. 

unlike the Brownian-coagulation mechanism; and thus 
it should be relevant to spinodal decomposition in quan¬ 
tum systems, which often takes place at extremely low 
temperatures leHii]. We note that in such a condition, 
quantum diffusion associated with tunnelling mass trans¬ 
fer can be operative but there is little thermal motion of 
droplets. Our study shows that the coupling between the 
composition and velocity fields alone can lead to domain 
coarsening even with little thermal noise. The relative 
importance of the velocity field to the composition field 
monotonically increases with an increase in the volume 
symmetry between the two phases; and thus the primary 
mechanism of coarsening changes from the LSW, to our 
mechanism, and to the Siggia’s mechanism with an in¬ 
crease of 4> from 0 to 1 / 2 . We hope that our mechanism 
will provide a novel perspective for spinodal decompo¬ 
sition of classical and quantum fluid systems in various 
fields. 


Methods 

Basic dynamic equations of model H. For our simu¬ 
lations, we employ the following kinetic equations known 
as a fluid model, or model H [TlElin]: 


dcf) 

dt 


- V V?!. + 


sm) 

5(j) 




( 2 ) 


p— + p(y ■ V)v = F 0 - Vp + + C, (3) 

where (j) is the composition, v is the fluid velocity, p is a 
part of pressure, p is the density, p is the viscosity, L is the 
kinetic coefficient, and [3 = 1/k^T (k^\ the Boltzmann 
constant). The pressure p is determined to satisfy the 
incompressiblity condition V • v = 0 . Since our primary 
concern is the dynamical coupling between the compo¬ 
sition and flow field, we do not care about the detailed 
properties of the mixture and thus adopt the following 
standard Ginzburg-Landau free energy functional: 

m = J +L|V<^|2], (4) 

where 7 = 7 o(T — Tc) and 70 , u and K are posi¬ 
tive constants. In Eq. is the thermodynamic 

force density acting on the fluid due to spatial inho¬ 
mogeneity of the composition field and = — 0 V/i = 
—V 7 r+/cBTLr(/)V^V 0 (tt: osmotic compressibility), where 
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/i = 51-L/5(j) is the chemical potential. Here 0 and C are 
the thermal noise terms in composition and force, which 
satisfy the following fluctuation-dissipation relation, re¬ 
spectively mill- 

(6>(r, t)0(Y , t )) = —2L\/^S{r — r )6{t — t (5) 

and 

(Ci(r,i)Ci(r = -2k-B.Tr]5ijV‘^5{v - r )8{t - t). (6) 

Here we rewrite Eqs. and by scaling the space 
and time unit by the correlation length of critical compo¬ 
sition fluctuations ^ and the characteristic 

lifetime of fluctuations /{\^\L). Then we define 

new variables, f = r/^ (r being the position vector), 
t = t/r^ and v = (^/r)v. The composition is normalized 
as 0 = <p/<pei where is the equilibrium 

composition of the coexisting phase. The scaled equa¬ 
tion which corresponds to Eq. (§ is then 


Numerical scheme to solve the above set of equa¬ 
tions properly with both composition and force 
noises and the details of simulations. Here we ex¬ 
plain how thermal composition and velocity fluctuations 
can be incorporated in our simulations. Although we 
have written Eqs. (7) and (8) as two separate equations, 
they can be unified into only one kinetic equation in the 
Stokes regime, since the velocity field is set solely by the 
composition field. Substituting Eq. (8) into Eq. (7), the 
dynamic equation of composition is given by [H [35] 


d(i>{r, i) 
di 



+ d^ir,i), 


( 12 ) 


where the kinetic coefficient is given as 


I/(f, f ) = V^(5(f - f ) + V0(f) • T(f - f ) • V0(f ) 


and the random noise term 0^{t) = — v^(f) • ^ 

satisfies the following fluctuation-dissipation relation 


at 

and Eq. (§ with the Stokes approximation becomes 
Vi = jdr + -'87^^)'87(i) + Cj, (8) 


where Tij is the Oseen tensor given by 
^ ~ Snr^ f + f3 


( 9 ) 


0 and ( are scaled thermal noises which satisfy the fol¬ 
lowing fluctuation-dissipation relations: 


{0{r,i)0{r ,i )) =—BV‘^6{r — r )6{i — i ), 


( 10 ) 


and 

{Ci{r,i)Cj{r ,i')) = -{B/A)5ijV‘^8{i-r)5{i-i'). (11) 

Here we note that (j/{k^T) = (2^/^/3)|7|^(/)g, which is the 
interface tension in an equilibrium two-liquid coexistence 
state. 


Simulations without thermal noise. We performed 
simulations without thermal noise to show that droplet 
coarsening takes place even without thermal noise (z.e., 
S = 0). We note that the Brownain coagulation mecha¬ 
nism fully relies on the presence of thermal force noise. 
Simulations without thermal noise are also useful to see 
the spatial distribution of physical quantities such as the 
chemical potential and the fluid velocity without suffer¬ 
ing from thermal noise. Eor these simulations, we in¬ 
troduce Gaussian random noise into 0 with an intensity 
5^ = 0.01 around the averaged composition 
(z.e., $ = 0.25) as an initial condition only at t = 0, 
switch off noise for t > 0, and follow the pattern evolu¬ 
tion process. 


(6>^(f,t)6>^(f )) = 2I/(f,f )S{i — t ). (13) 

Here the random velocity noise v^ is given by v^(f, t) = 
f dr T(f — f )C(f)- Since the kinetic coefficient of the 
above equation depends upon the composition field 0, 
the noise term is multiplicative and should be treated via 
the Stratonovich interpretation [53|. Thus we adopted a 
Stratonovich scheme with semi-implicit terms developed 
by Camley and Brown [44] under a periodic boundary 
condition with a staggered grid. This scheme was orig¬ 
inally conceived to study the phase separation process 
in a two-dimensional membrane. We modify their ap¬ 
proach by using the Oseen tensor appropriate for a three- 
dimensional (3D) Newtonian fluid in the Stokes regime. 
Here it is worth noting that 3D simulations of model H 
with thermal noises were also performed by a Lattice- 
Boltzmann method [46]. The system size was set to 
256^ lattices, which is large enough to study the late 
stage of phase separation. We set the grid size to be 
Ax = Ay = Az = 1. The volume fraction of the minor¬ 
ity phase 4> is set to 0.25 and the composition is set to 
spatially homogeneous at t = 0. 

Simulations of the motion of a pair of droplets. We 

place two droplets whose interfacial profile is described 
by a step function at t = 0 and then follow the process 
of equilibration by solving the model H equations. In 
this simulation, the fluidity parameter A is 50 and the 
thermal fluctuations are not included (ie., S = 0). Here 
the large fluidity is used to access a large separation be¬ 
tween the interfacial thickness and the droplet size within 
a limited simulation time as well as to see hydrodynam¬ 
ics effects clearly. Note that the process is much slower 
without thermal noises. The simulation box size is 192^. 

Analysis. We analyse the sizes and positions of droplets 
by counting the number of lattice points belonging to 
each droplet and calculating the positions of the centre 
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of mass of the lattice points after binarization of the com¬ 
position field. We also obtain the velocity of a droplet by 
calculating the average of velocity over lattice points in¬ 
side the droplet. The structure factor S{k) is obtained by 


spherically averaging the power spectrum of the Fourier 
transformation of 0(f). The characteristic wavenumber 
ki is calculated as = f dk kS{k)/ f dk S{k). 
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